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I. INTRODUCTION 



The Maxwell's equations in vacuum can be written as 
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where E is the electric field vector, B is the magnetic 
induction vector, p e is the density field of electric charges, 
£q is the dielectric constant of vacuum, /xq is magnetic 
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permeability of vacuum, t is time, V = + j 
is the Hamilton operator. 

The main purpose of this paper is to derive the above 
mentioned Maxwell's equations in vacuum based on a 
continuum mechanics model of vacuum and a singularity 
model of electric charges. 

The motive of this paper is to seek a mechanism of 
electromagnetic phenomena. The reasons why new me- 
chanical models of electromagnetic field are interesting 
may be summarized as follows. 

Firstly, there exists some electromagnetic phenomena 
that could not be interpreted by the present theories of 
electromagnetic field, e.g.,the spin of electrons [JQ, the 
Aharonov-Bohm effect [3|, |4J. New theories of of electro- 
magnetic phenomena may view these problems from new 
angles. 

Secondly, there exists some inconsistencies and inner 
difficulties in the classical electrodynamics, e^g., the in- 
adequacy of Lieenard-Wiechert potentials [a Is ■ New 
theories of electromagnetic phenomena may overcome 
such difficulties. 

Thirdly, there exists some divergence problems in 
quantum electrodynamics [§]. In Dirac's words, 'I must 
say that I am very dissatisfied with the situation, be- 
cause this so-called good theory does involve neglecting 
infinities which appear in its equations, neglecting them 



in an arbitrary way. This is just not sensible mathemat- 
ics.' New theories of electromagnetic phenomena may 
open new ways to solve such problems. 

Fourthly, since quantum theory shows that vacuum 
is not empty and has physical effects, e.g., the Casimir 
effect H US EH, El, it is valuable to reexamine the old 
concept of electromagnetic aether. 

Fifthly, from the point view of reductionism, Maxwell's 
theory of electromagnetic field can only be regarded as 
a phenomenological theory. Although Maxwell's theory 
is a field theory the concept of field is different from 
that of continuum mechanics ID, 14 , 15, li[ because of 
the absence of a continuum. Thus, from the point view 
of reductionism, the mechanism of electromagnet ic p he- 
nomena remains an unsolved problem in physics [171 ] . 

Sixthly, one of the puzzles in physics is the problem 
of dark matter and dark energy (refer to, for instance, 



El, 0, HE HI Hi Hi Hi)- New theories of 
electromagnetic phenomena may provide new ideals to 
attack this problem. 

Finally, one of the tasks of physics is the unification of 
the four fundamental interactions in the universe. New 
theories of electromagnetic phenomena may shed some 
light on this puzzle. 

To conclude, it seems that new considerations on elec- 
tromagnetic phenomena is needed. It is worthy keeping 
an open mind with respect to all the theories of elec- 
tromagnetic phenomena before the above problems been 
solved. 

Now let us briefly review the long history of the me- 
chanical interpretations of electromagnetic phenomena. 

According to E. T. Whittaker[l7(, Descartes was the 
first to bring the concept of the aether into science by sug- 
gesting that it has mechanical properties according to E. 
T. Whittaker 17] . Descartes believed that every physi- 
cal phenomenon could be interpreted in the construction 
of a mechanical model of the universe. William Watson 
and Benjamin Franklin introduced the one- fluid theory of 
electricity independently in 1746 17]. Henry Cavendish 
attempted to explain some of the principal phenomena of 
electricity by means of an elastic fluid in 1771 [17|. Not 
contented with the above mentioned one-fluid theory of 
electricity, du Fay, Robert Symmer and C. A. Coulomb 



1 



developed a two-fluid theory of electricity from 1733 to 
1789 0. 

Before the unification of electromagnetic phenomena 
and light phenomena by Maxwell in 1860s, light phenom- 
ena were also studied independently based on Descartes' 
scientific research program of the mechanical theory of 
nature. John Bernoulli introduced a fluidic aether the- 



abandoned the concept of aether 27J. However, Ein- 
stein's assertion did not cease the exploration of aether 
(refer to for instance, [13, 0, 0, S EE (H [H, [34, 



ory of light in 1752 17]. Euler believed that all electrical 
phenomena is caused by the same aether that propagates 
light. Furthermore, Euler attempted to explain gravity 
in terms of his single fluidic aether [l7j|. 

In 1821, in order to explain polarisation of light, A. 
J. Frensnel proposed an aether model which is able to 
transmit transverse waves. After the advent of this suc- 
cessful transverse wave theory of light of Frensnel, those 
imponderable fluid theories were abandoned. Frensnel's 
dynamical theory of a luminiferous aether had an impor- 
tant influence on the mechanical theories of nature in the 
nineteenth century [ItJ ■ Inspired by Frensnel's luminifer- 
ous aether theory, numerous dynamical theories of elastic 
solid aether were established by Stokes, Cauchy, Green, 
MacCullagh, Boussine sq . Riemann and William Thom- 
son, see, for instance, [171 ]. 

Thomson's analogies between electrical phenomena 
and elasticity helped to inspire James Clark Maxwell 
established a mechanical model of electrical phenomena 
[l7j ]. Strongly impressed by Faraday's theory of lines of 
forces, Maxwell compared the Faraday's lines of forces 
with the lines of flow of a fluid. In 1861, in order to ob- 
tain a mechanical interpretation of electromagnetic phe- 
nomena, Maxwell established a mechanical model of a 
magneto-electric medium. Maxwell's magneto-electric 
medium is a cellular aether, looks like a honeycomb. 
Each cell of the aether consists of a molecular vortex 
surrounded by a layer of idle- wheel particles. In a re- 
markable paper published in 1864, Maxwell established 
a group of equations, which were named after his name 
later, to describe the electromagnetic phenomena. 

In 1878, G. F. FitzGerald compared the magnetic force 
and the velocity in a quasi-elastic solid of the type first 
suggested by MacCullagh 17j. FitzGerald's mechanical 
model of the electromagnetic aether were studied by A. 
Sommerfeld, by R. Reiff and by Sir J. Larmor later [171 ]. 

Because of some dissatisfactions with the mechanical 
models of the electromagnetic aether and the success 
of the theory of electromagnetic field, the mechanical 
world view was replaced by the electromagnetic world 
view gradually. Therefore, the concepts of a luminif- 
erous aether and an elastic solid aether were replaced 
by the concepts of a electromagnetic aether or the elec- 
tromagnetic field. This paradigm shift in scientific re- 
search was attributed to many scientists, including Fara- 
day, Maxwell, Sir J. Larmor, H. A. Lorentz, J. J. Thom- 
son, H. R. Hertz, Ludwig Lorenz, Emil Wiechert, Paul 
Drude, Wilhelm Wien, etc., see, for instance, (l7| . 

In a remarkable paper published in 1905, Einstein 



351 . 1361 . |37j). Einstein changed his attitude later and 
introduced his new concept of ether 38, [39[. In 1979, 
A. A. Golebiewska-Lasta observed the similarity between 
the electromagnetic field and the linear dislocation field 
[jjj ]. V. P. Dmitriyev have studied the similarity be- 
tween the electromagnetism and the linear elasticity since 
1992[H[3i,[33,[3]. H. Marmanis established a new the- 
ory of turbulence based on the analogy between electro- 
magnetism and turbulent hydrodynamics in 1998 34j. D. 
J. Larson derived the Maxwell's equations from a simple 
two-component solid-mechanical aether in 1998 [33j. D. 
Zareski gave an elastic interpretation of electrodynamics 
in 2001 [36|. I regret to admit that it is impossible for me 
to mention all the works related to this field in history. 

Inspired by the above mentioned works, we show that 
the Maxwell's equations of electromagnetic field can be 
derived based on a continuum mechanics model of vac- 
uum and a singularity model of electric charges. 



II. CLUES OBTAINED FROM DIMENSIONAL 
ANALYSIS 

According to Descartes' scientific research program 
which was based on his mechanical view of nature, the 
electromagnetic phenomenon must be and can be inter- 
preted based on the mechanical motions of aether parti- 
cles. 

Therefore, all the physical quantities appearing in the 
theory of electromagnetic field must be mechanical quan- 
tities. 

Thus, in order to establish a successful mechanical 
model of electromagnetic field, we should undertake a 
careful dimensional analysis (refer to, for instance, [ilj]). 
of electric field vector E, magnetic induction vector B, 
the density field of electric charges p q , the dielectric con- 
stant of vacuum eo, magnetic permeability of vacuum /io, 
etc. 

It is known that the Maxwell's equations (Q][I]) in vac- 
uum can also be expressed as [l[ 
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where cf> is the scalar electromagnetic potential, A is the 
vector electromagnetic potential, V = + j^p + 
is the Laplace operator. 

Noticing the similarity between the Eq.© and the 
equation ([38|1 of momentum conservation of elastic solids, 
it is natural to speculate that the dimension of the vec- 
tor electromagnetic potential A of aether has the same of 
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the dimension of the displacement vector u of an elastic 
solid. 

In 1846, W. Thomson compared electric phenomena 
with elasticity. He pointed out that the elastic displace- 
ment u of an incompressible elastic solid is an possible 
analogy with the vector electromagnetic potential A |17|. 

Inspired by this clue, let us set out to investigate in 
this direction further in the following sections. 



III. A VISCO-ELASTIC CONTINUUM MODEL 
OF VACUUM 

The purpose of this section is to establish a visco- 
elastic continuum mechanical model of vacuum. 

In 1845-1862, Stokes suggested the aether might be- 
have like a glue- water jelly |42l. 0, He believed that 
the aether would act like a fluid with respect to transla- 
tional motion of large bodies through it, but would still 
posses elasticity to produce small transverse vibrations. 

Following Stokes, we propose a visco-elastic continuum 
model of vacuum. 

Hypothesis 1 Suppose the universe is filled with a kind 
of continuously distributed material. 

This material may be called aether for convenience. 

In order to establish a continuum mechanical theory of 
aether, we must introduce some hypotheses based on ex- 
perimental data of the macroscopic behavior of vacuum. 

Hypothesis 2 We suppose that the all the mechanical 
quantities of the aether under consideration, such as den- 
sity, displacements, strains, stresses, etc., are piecewise 
continuous functions of space and time. Furthermore, 
we suppose that the material points of the aether remain 
one-to-one correspondence with the material points before 
deformation happens. 

Hypothesis 3 We suppose that the material of the 
aether under consideration is homogeneous, that is ^ = 
= ff = af = Oj where p is the density of aether. 

Hypothesis 4 Suppose that the deformation processes 
of aether are isothermal processes. We neglect the ther- 
mal effects. 

Hypothesis 5 Suppose the deformation processes is not 
influenced by the gradient of stress tensor. 

Hypothesis 6 We suppose that the material of the 
aether under consideration is isotropic. 

Hypothesis 7 We suppose that the deformaton of the 
aether under consideration is small. 

Hypothesis 8 We suppose that there are no initial 
stress and strain in the body under consideration. 



When aether is subjected to a set of external forces, 
the relative positions of the aether particles forming the 
body changes. 

In order to described the deformation of the aether, let 
us introduce a Cartesian coordinate system {o, x, y, z} or 
{o, x±, X2, £3} which is static relative to the aether. Now 
we may introduce the definition of displacement vector u 
of every point in the aether as 



ro, 
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where ro is the position of the point before the deforma- 
tion, r is the position of the point after the deformation. 

The displacement vector may be written as u = u\i + 
U23 + 113k. or u = ui + uj + u>k, where i, j, k are the three 
unit vectors along the three coordinates. 

The gradient of the displacement vector u is the rela- 
tive displacement tensor Uij — 

We can decompose the tensor Uij into two parts, sym- 
metric eg and skew-symmetric u>ij (refer to, for instance, 

USUI)- 
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-(Uij + U j}i ) + -{Ui,j - Uj.i) = Si 
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Hj = -z(uij -Uj,i). (9) 



The symmetric tensor Ey represents pure deformation 
of the body at a point and is called strain tensor (refer 
to, for instance, 14, 45, 46[). The matrix form and the 



indicial notation of strain tensor Sij are 
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The strain-displacements equations can be obtained 
from Ea.([TO|) 
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For convenience we can introduce the definitions of 



mean strain e m and strain deviator en as 
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When the aether is deformed, internal forces arise due 
to the deformation. The the indicial notation of the stress 
tensor er^ is 



Oil 012 013 
021 022 023 
V0"31 032 033/ 



(14) 



For convenience we can introduce the definitions of 
mean stress a Ti 



and stress deviator Sij as 
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Now let us turn to study the constitutive relation. 

An elastic Hooke solid responds instantaneously with 
respect to a external stress. A Newtonian viscous fluid 
responds to a shear stress by a steady flow process. 

In 19th century, people began to notice that some ma- 
terials showed time dependence in their elastic response 
with respect to external stress. When materials like 
pitch, gum rubber, polymeric materials, hardened ce- 
ment and even glass were loaded, an instantaneous elas- 
tic deformation was followed by a continuous slow flow 
or creep. 

Now this time-dependent response is called viscoelas- 
ticity (refer to, for instance, [13, S, Mate- 
rials exhibits both an instantaneous elastic elasticity 
and creep characteristics is called viscoelastic materials 
HI, . Viscoelastic materials were studied long ago by 
Maxwellfii, [H, [H3], Kelvin, Voigt, BoltzamannJIl 0, 
51j . etc. 

Inspired by these contributors, we propose a visco- 
elastic constitutive relation of the aether. 

It is natural to speculate that the constitutive relation 
of the aether may be a combination of the constitutive 
relations of the Hooke-solid and the Newtonian-fluid. 

For the Hooke-solid, we have the follo wing the general- 
ized Hooke's law (refer to, for instance, [l4 liBl liil IHij] ) . 
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where Sij is the Kronecker symbol, a m is the mean stress, 
where Y is the Yang modulus, v is the Poisson ratio, 
G is the shear modulus, A is Lame constant, 9 is the 
volume change coefficient. The definition of 9 is 9 = 
£11 + £22 + £33 - + ^ + -57- 

The generalized Hooke's law Eq. (|T7|) can also be writ- 
ten as 15211 
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where Sjj is the stress deviator, 



is the strain deviator. 



For the Newtonian-fluid, we have the following consti- 
tutive relation 



de 
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where Sij is the stress deviator, deij/dt is the strain rate 
deviator, r\ is the dynamic viscocity. 

Since the aether behaves like the Hooke-solid during 
very short time intervals, we therefore differentiate both 
sides of Eq. (fTB|) and obtain 
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A combination of Eq. ((20)) and Eq. (fT9|) gives 
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The materials which behave like Eq. pTj) are called 
Maxwell-liquid since Maxwell established such a consti- 
tutive relation in 1868 (refer to, for instance, 47, 48|, 4^, 

Eq. (|2~lj) is valid only in the case of infinitesimal defor- 
mation because of the presence of the derivative with 
respect to time. Oldroyd recognized that we need a 
special definition of derivative operation in order to 
satisfy the principle of material frame indifference or 
objectivity 0,[H3|. Unfortunately, there is no unique def- 
inition of derivative operation fulfil the the principle of 
objectivity presently 



As an enlightening example, let us recall the descrip- 
tion in [47J about a simple shear experiment. We suppose 
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where cr t is the shear stress, ej is the shear strain. 
Therefore, Eq. (f2"Tj) becomes 
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Integration of Eq. |23|) gives 
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If the shear deformation is kept constant, i.e. de t /dt 
0, we have 



a t = a e 1 



(25) 



Eq. p5| shows that the shear stresses remain in the 
Maxwell-liquid and are damped in the course of time. 

We see that 77/G must have the dimension of 
time. Now let us introduce the following definition of 
Maxwellian relaxation time r 
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G' 
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Therefore, using Eq.([26|). Eq. (|2Tj) becomes 



Noticing Eq.([28]), Eqs. (|30H32[) can also be written as 



dsij 
~dt 



Hp- 
= 2G^ 
dt 



(27) 



Now let us introduce the following hypothesis 



Hypothesis 9 Suppose the constitutive relation of the 
aether satisfies Eg. H21\ ). 



Now we can derive the the equation of momentum con- 
servation based on the above hypotheses [9] 

Let T be the characteristic time scale of a observer 
of the aether. When the observer's time scale T is the 
same order of the period of wave motion of light, the 
Maxwellian relaxation time r is a relatively a large num- 
ber. Thus, the first term of Eq. ([27|) may be neglected. 
Therefore, the observer concludes that the strain and the 
stress of the aether satisfies the generalized Hooke's law. 

The generalized Hooke's law (| 1 T|) can also be written 
as [illlsl 
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where A = Yv/[(1 + v)(\ — 2v)\ is Lame constant, 6 
is the volume change coefficient. The definition of 9 is 
V - en + £22 + £33 - -57 + + -gj- 
The following relationship are useful 
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where K is the volume modulus. 

It is known that the equations of momentum conser- 
vation are (refer to, for instance, 14, 45, 4(| H2, 54. HB|). 



(9(7n 


(9(712 


9(713 


Yf x 


d 2 u 


dx 


dy 


dz 


= P W 


9(721 


9(7 22 


9(723 

V 


Yfy 


d 2 v 


dx 


dy 


dz 


= P W 


do-31 


da 32 


9(733 


Yh 


d 2 w 


dx 


dy 


dz 


= P dt 2 



(30) 

(31) 
(32) 



where f x , f y and f z are the volume force density exerted 
on the aether. 

The tensor form of the equations (|30ti32|) of momentum 
conservation can be written as 
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Using Eq. pTjl . Eqs. (|34j|35|) can also be expressed by 
means of displacement u 

9 / du dv dw\ 



GV^ + (G + A)— l— + — + —)+f x = p—, 
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(37) 

The vector form of the above equations (1371) can be 



inc vector iorm 01 tnc above equations \\d i\) can be 
written as (refer to, for instance, 1J, 45, 46|, 521.1541 IHBT]). 
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When there are no body force in the aether, Eqs. 
reduces to 



GV 2 u+(G + A)V(V-u) 



9 2 u 



(39) 



From a theorem of Long 4^, 5(|, there exist a scalar 
function <f> and a vector function R such that u is repre- 
sented by 

u = \/tp + V x R (40) 
and <fi and R satisfy the following wave equations 
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where q is the velocity of longitudinal waves, Ct is the 
velocity of transverse waves. The definitions of these two 
elastic wave velocities are (refer to, for instance, 45, 4(1 



5J,l55j) 
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ip and R is usually called the scalar displacement 
potential and the vector displacement potential respec- 
tively. 
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IV. DEFINITION OF POINT SOURCE AND 
SINK 

If there exists a velocity field which is continuous and 
finite at all points of the space, with the exception of 
individual isolated points, then these isolated points are 
called velocity singularities usually. Point sources and 
sinks are examples of velocity singularities. 

Definition 10 Suppose there exists a singularity at point 
Pq = (xq, yoi Zq) in a continuum. If the velocity field of 
the singularity at point P = (x,y,z) is 



v(x,y,z,t) 



Q } 

Anr 2 



(44) 



where r = y/(x — xo) 2 + (y — yo) 2 + (z — zo) 2 , f denotes 
the unit vector directed outward along the line from the 
singularity to the point P = (x,y,z), then we call this 
singularity a point source if Q > or a point sink if 
Q < 0. Q is called the strength of the source or sink. 

Suppose a static point source with strength Q locates 
at the origin (0,0,0). In order to calculate the volume 
leaving the source per unit time, we may enclose the 
source with an arbitrary spherical surface S with radius 
a. A calculation shows that 



u • ndS — 



Q 



w r • ndS = Q, 



(45) 



where n denotes the unit vector directed outward along 
the line from the origin of the coordinates to the field 
point (a;, y, z). Equation (j45]) shows that the strength Q of 
a source or sink evaluates the volume of the fluid leaving 
or entering a control surface per unit time. 

For the case of continuously distributed point sources 
or sinks, it is useful to introduce the following definition 
of volume density p s of point sources or sinks, 



= lim , 



(46) 



where Ay is a small volume, AQ is the sum of strengthes 
of all the point sources or sinks in the volume AV\ 



V. A POINT SOURCE AND SINK MODEL OF 
ELECTRIC CHARGES 

The purpose of this section is to propose a point source 
and sink model of electric charges. 

Let T be the characteristic time of a observer of a elec- 
tric charge in the aether. We may suppose that the 
observer's time scale T is very large compares to the 
Maxwellian relaxation time r. So the Maxwellian relax- 
ation time r is a relatively a small number and the stress 
deviator Sij changes very slowly. Thus, the second term 



in the left of Eq. (j27j) may be neglected. According to this 
observer, the constitutive relation of the aether may be 
written as 



(47) 



Therefore, the observer concludes that the aether be- 
haves like the Newtonian-fluid under his time scale. 

In order to compare fluid motion with electric fields, 
Maxwell introduced an analogy between sources or sinks 
and electric charges (l7j ]. Inspired by Maxwell, we may 
introduce the following 

Hypothesis 11 Suppose all the electric charges in the 
universe are sources or sinks in the aether. We define a 
source as a negative electric charge. We define a sink as 
a positive electric charge. The electric charge quantity q e 
of a electric charge is defined as 



-k Q pQ, 



(48) 



where p is the density of the aether, Q is called the 
strength of the source or sink, kg is a positive dimen- 
sionless constant. 



A calculation shows that the mass 
charge is changing with time as 



dm 
~dt 



-PQ 



Qe_ 



of a electric 



(49) 



where q e is the electric charge quantity of the electric 
charge. 

We may introduce a hypothesis that the masses of elec- 
tric charges are changing so slowly relative to the time 
scaler of human beings that they can be treated as con- 
stants approximately. 

For the case of continuously distributed electric 
charges, it is useful to introduce the following definition 
of volume density p e of electric charges 



Pp = lim — — 

y av^o AV 



(50) 



where AV is a small volume, Aq e is the sum of strengthes 
of all the electric charges in the volume AV. 
From Eq.([4"6"l). Eq.(T48]) and Eq.pj]). we have 

Pe = -kqpps (51) 



VI. DERIVATION OF THE MAXWELL'S 
EQUATIONS IN VACUUM 

The purpose of this section is to derive the Maxwell's 
equations based on the above visco-elastic continuum 
model of vacuum and the singularity model of electric 
charges. 
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Now let us to derive the continuity equation of the 
aether from mass conservation. Consider an arbitrary 
volume V bounded by a closed surface S fixed in 
space. Suppose there are electric charges continuously 
distributed in the volume V. The total mass in volume 
V is 



M 



pdV, 



(52) 



where p is the density of the aether. 

The rate of increase of the total mass in volume V is 



ME- 



(53) 



Using the Ostrogradsky-Gauss theorem (refer to, for 
instance, [IE M, S M, H ED, H3), the rate of mass 
outflow through the surface S is 



p(u ■ n)dS 



V • (pu)dV, 



where v is the velocity field of the aether. 
The definition of the velocity field v is 



(54) 
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Using Eq. (|5Tj) , the rate of mass created inside the vol- 
ume V is 



PPsdV 



kQ 



-dV. 



(56) 



Now according to the principle of mass conservation, 
and making use of Eq. (f53")) . Eq. (f54")l and Eq. (j56"|) , we have 



V • (pv)dV (57) 



Since the volume V is arbitrary, from Eq. (|57p we have 



dp 
dt 



V ■ (pv) 



£e_ 

kQ 



(58) 



According to Hypothesis^! the aether is homogeneous, 
that is|£ = §£ = !f = ff=0. Thus, Eq.® becomes 



Pe 

kQp' 



(59) 



According to Hypothesis [TT1 and Eq. (|4"9")l . the masses of 
positive electric charges are changing since the strength 
of a sink evaluates the volume of the aether entering the 
sink per unit time. Thus, the momentum of a volume 
element AV of the aether containing continuously dis- 
tributed electric charges moving with an average speed 
v e is changing. The increased momentum AP of the 
volume element AV during a time interval At is the de- 
creased momentum of continuously distributed electric 



charges contained in the volume element AV during a 
time interval At, that is, 



AP = p(p s AVAt)v e = --p-AVA£v e (60) 

kQ 

Therefore, the equation of momentum conservation 
Eq. ([38]) of the aether should be changed as 



GV 2 u + (G + A)V(V ■ u) + f = p 



d 2 u _ p e v e 
dt 2 k 



Q 



(61) 



In order to simplify the Eq. (|6ip . we may introduce an 
additional hypothesis as 

Hypothesis 12 We suppose that the aether is incom- 
pressible, that is 8 = 0. 

From Hypothesis [TSJ we have 

du dv dw 
V-u= — + — + — = 6 = 62 

dx dy dz 

Therefore, the vector form of the equation of momen- 
tum conservation Eq. (|61|) reduces to the following form 



d 2 u p e v e 



GV^u + f=p-- T - 



dt 2 k Q ■ 



(63) 



According to the Stokes-Helmholtz resolution theorem 
(refer to, for instance, 0, IH^I). which states that ev- 
ery sufficiently smooth vector field may be decomposed 
into irrotational and solenoidal parts, there exist a scalar 
function <j> and a vector function <j) and a R such that u 
is represented by 



u = Vt/> + V x R. 



Now let us introduce the definitions 



(64) 



V<£ = fc E — (WO, A = fc B VxR, (65) 



du 

V = -k E — , B = fc B Vxu, (66) 

where <p is the scalar electromagnetic potential, A is the 
vector electromagnetic potential, E is the electric field 
intensity, B is the magnetic induction, fcg is a positive 
dimensionless constant. 

From Eq.ljBip. Eq.^ and Eq.®, we have 



dA 

E = -V(f>~ — , B=VxA 
dt 



and 



V x E = 
V-B = 0. 



dB 

~dt' 



(67) 



(68) 
(69) 
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Based on Eg. ([65)1 and noticing 

V 2 u = V(V • u) - V x (V x u), 
V 2 A = V(V • A) - V x (V x A), 

and V • u = 0, V • A = 0, we have 



fc K V 2 u = V 2 A. 



Therefore, using Eg. ([72]) . Eq. (|63|) becomes 
G 2 a , c d 2 u p e w e 



Using Ea.([7T]). Eq.([75|) becomes 

G _ ,_ . , _ d 2 u p e v e 
- — V x V x A + f = p— - 

k E dt 2 k Q 



Now using Eq. (p7|) , Eg. (1741 becomes 

G „ „ p 9E p e v e 
V x B + f = — t-—^- 

k E k E dt k Q 



(70) 
(71) 



(72) 



(73) 



(74) 



(75) 



It is natural to speculate that there are no other body 
forces or surface forces exerted on the aether. Thus, we 
have f = 0. Therefore, Eq. ([75|) becomes 



k Q G 



V x B 



k Q pdE 



k E k E dt 

Now let us introduce the following definitions 
kQp 1 _ k Q G 



eo 



k E 



Mo 



k E 



(76) 



(77) 



Therefore, Eq. ([76|) becomes 
1 



„ . <9E 

— V x B = j + eo-^r 
Po dt 



(78) 

Noticing Eq.® and Eq.([77|l. Eq.(|59]l becomes 

V-E=— . (79) 



Pe 

eo 



Now we see that Eq.([55]). Eq, 1(55)1. Eq.® and Eq. 1(751 
coincide with the Maxwell's equations (PHJ). 



VII. MECHANICAL INTERPRETATION OF 
THE ELECTROMAGNETIC WAVE 

It is known that from the Maxwell's equations (DID), we 
can obtain the following equations (refer to, for instance, 



1 <9 2 E 1 dj 
V 2 E — = -V Pe + 

p e dt 2 e dt 

„ 2tt 1 5 2 H 1 „ . 

V H -TJT = Vx J- 

/ioeo po 



(80) 
(81) 



Eg. ((80)1 and Eg. ((81)) are the electromagnetic wave 
equations with sources in the aether. In the source free 
region where p e = and j = 0, this equations reduce to 
the following equations 



V 2 E- 



1 <9 2 E 



p e dt 2 

V 2 H-^^ 

p e dt 2 



= 0, 



= 0. 



(82) 
(83) 



Eg. ((82)) and Eg. (|83|) are the electromagnetic wave equa- 
tions without sources in the aether. 

From Eq. ((82l . Eq.([g5j> and Eq.([77j). we see that the 
velocity Co of electromagnetic waves in vacuum is 



1 



co = 



V^oeo 



(84) 



From Eg. (f43)) and Eg . (|84|) . we see that the velocity Cq 
of electromagnetic waves in the vacuum is the same as the 
velocity ct of the transverse elastic waves in the aether. 

Now we may regard electromagnetic waves in the vac- 
uum as the transverse waves in the aether. This idea was 



first introduced by Fresnel in 1821 [17|. 



VIII. CONCLUSION 

It is an old idea that the universe may be filled with 
a kind of continuously distributed material which may 
be called aether. Following Stokes, we propose a visco- 
elastic constitutive relation of the aether. Following 
Maxwell, we propose a fluidic source and sink model 
model of electric charges. Thus, the Maxwell's equations 
in vacuum are derived by methods of continuum mechan- 
ics based on this continuum mechanical model of vacuum 
and the singularity model of electric charges. 



IX. DISCUSSION 

There exists some interesting theoretical, experimen- 
tal and applied problems in the fields of continuum me- 
chanics, the classical electrodynamics, quantum electro- 
dynamics and other related fields involving this theory of 
electromagnetic phenomena. It is an interesting task to 
generalize this theory of electromagnetic phenomena in 
the static aether to the case of electromagnetic phenom- 
ena of moving bodies. 
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